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ABSTRACT 

We present a method for selecting high-redshift type la supernovae (SNe la) 
located via rolling SN searches. The technique, using both color and magnitude 
information of events from only 2-3 epochs of multi-band real-time photometry, 
is able to discriminate between SNe la and core collapse SNe. Furthermore, for 
the SNe la, the method accurately predicts the redshift, phase and light-curve 
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parameterization of these events based only on pre-maximum-light data. We 
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demonstrate the effectiveness of the technique on a simulated survey of SNe la 
and core-collapse SNe, where the selection method effectively rejects most core- 
collapse SNe while retaining SNe la. We also apply the selection code to real-time 
data acquired as part of the Canada-France-Hawaii Telescope Supernova Legacy 
Survey (SNLS). During the period May 2004 to January 2005 in the SNLS, 
440 SN candidates were discovered of which 70 were confirmed spectroscopically 
as SNe la and 15 as core-collapse events. For this test dataset, the selection 
technique correctly identifies 100% of the identified SNe II as non-SNe la with 
only a 1-2% false rejection rate. The predicted parameterization of the SNe la 
has a precision of |Az| /(l + -2 sp cc) < 0.09 in redshift, and ±2-3 rest-frame days 
in phase, providing invaluable information for planning spectroscopic follow-up 
observations. We also investigate any bias introduced by this selection method 
on the ability of surveys such as SNLS to measure cosmological parameters (e.g., 
w and Om), and find any effect to be negligible. 

Subject headings: surveys - supernovae: general - distance scale 

1. Introduction 

Recent surveys for high-redshift type la supernovae (SNe la) have established their 
utility as calibrated standard candles suitable for use as cosmological distance indicators. The 
"first generation" of SN la cosmology programs (Perlmutter et al. 1997, 1998; Garnavich et al. 
1998; Riess et al. 1998; Schmidt et al. 1998; Perlmutter et al. 1999) developed a systematic 
approach to high-redshift SN detection and analysis that led to Hubble diagrams ruling out a 
flat, matter-dominated Universe, revealing the presence of an unaccounted-for "dark energy" 
driving cosmic acceleration. 

One of the most pressing questions in cosmology now is: "What is the nature of the 
dark energy?" . There is a fundamental difference between a Cosmological Constant and other 
proposed forms of dark energy, and the distinction can be addressed by measuring the dark 
energy's average equation-of-state, w = p/p, where w — — 1 corresponds to a Cosmological 
Constant. Current measurements of this parameter (e.g. Knop et al. 2003; Tonry et al. 2003; 
Riess et al. 2004a) are consistent with a very wide range of dark energy theories (see Peebles 
& Ratra 2003, for a comprehensive review). The importance of improving measurements to 
the point where w — — 1 could be confirmed or excluded has led to a "second-generation" of 
supernova cosmology studies: large multi-year, multi-observatory programs benefiting from 
major commitments of dedicated time, e.g. the Supernova Legacy Survey at the Canada- 
France-Hawaii Telescope (Astier et al. 2005), and the ESSENCE project at the Cerro Tololo 
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Inter-American Observatory (Krisciunas et al. 2005). These "rolling searches" find and fol- 
low all varieties of SNe over many consecutive months of repeated wide-field imaging, with 
redshifts and SN type classification from coordinated spectroscopy. The success of these 
large-scale rolling searches directly depends on the degree to which SNe la can be selected 
at early phases in their light-curve from many other types of varying sources for further 
spectroscopic follow-up, and hence confirmation of their nature. With large numbers of 
candidates being found during each lunation, and a finite amount of 8m-class telescope spec- 
troscopic follow-up available, these spectroscopy-limited surveys must have reliable methods 
for extracting likely SN la events from the hundreds of other varying sources discovered. 

Previous studies have investigated methods by which SN can be identified from light- 
curve or host-galaxy photometric properties alone. Poznanski et al. (2002) used color-color 
diagrams and SN spectral templates at various redshifts to delineate regions of color-color 
space in which SNe of different types are usually located; a particular adaptation of this tech- 
nique is to identify SNe Ib/c (Gal- Yam et al. 2004). Dahlen & Goobar (2002) demonstrated 
that photometric properties of the host galaxies and SNe can be used to select high-redshift 
candidates in well-defined redshift intervals for spectroscopic follow-up. Riess et al. (2004b) 
used color cuts to differentiate SNe la at z > 1 from lower-redshift SNe II, while Barris et al. 
(2004) used photometric techniques to attempt to identify high-redshift SN candidates (with 
complete light-curves) as SNe la based on a spectroscopic redshift. Barris & Tonry (2004) 
introduced a redshift marginalization technique which allowed SNe la to be used without 
knowledge of a spectroscopic redshift with an impressively small increase in scatter about 
the z < 0.13 Hubble diagram; however this technique used full light-curves and an already 
obtained positive identification of the SN type. 

Many (though not all) of these techniques are designed to identify SNe either after a 
host galaxy spectroscopic redshift has been determined, or after a full light-curve has been 
obtained. The problem facing rolling-searches is different, namely differentiating SNe la 
from all other variable sources (e.g., SNe II, AGN, variable stars) based on only 2-3 epochs 
of real-time (i.e. perhaps imperfect) photometry, and predicting their redshifts, phases, and 
temporal evolution in magnitude. Redshift estimates allow spectrograph setups to be opti- 
mized (in terms of grism or grating central wavelength choice) for the candidate under study, 
to ensure that the characteristic spectral features of SNe la are always within the observed 
wavelength range to allow unambiguous identification. Phase predictions allow spectroscopic 
follow-up to be optimized by always observing SN candidates within a few days of maximum 
light, when the SNe have the greatest contrast over their host galaxies, and when spectral 
features in different types of SNe are most easily distinguished. Predicting the magnitude of 
a SN over the course of its light-curve is invaluable in ensuring that spectroscopic observa- 
tions are of a sufficient length to obtain a signal-to- noise (S/N) that allows an identification 
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to be obtained. 

In this paper, we present a new technique which is designed to select and predict the 
redshift, phase and light-curve parameterization of high-redshift SNe la. We show how it 
has been applied to real-time data from the Supernova Legacy Survey (SNLS) and used to 
guide and inform follow-up decisions. A plan of the paper follows. §2 gives a brief overview 
of SN la properties and describes the technique used for our photometric selection; §3 shows 
the results of our tests of the technique on synthetic datasets of SNe la and core collapse 
SNe. The technique is applied to candidates from the SNLS in § 4, the results of which are 
presented, compared with simulations, and discussed in §5. We summarise in §6. Throughout 
this paper all magnitudes are given in the AB photometric system (Oke & Gunn 1983), and 
we assume a cosmology of VL M = 0.30, Q,\ = 0.70, and w — — 1 where relevant. 

2. Photometric selection technique 

Given flux information about a candidate SN, we require a method that provides a light- 
curve fit and parameterization of the SN (assuming the candidate is a SN la), an estimate 
of the SN redshift, and some quality factor providing a measure of the likelihood that the 
candidate under study is indeed a SN la. 

2.1. Parameterizing Type la supernovae light-curves 

Raw observed peak magnitudes of spectroscopically normal SNe la have an r.m.s. dis- 
persion of ~ 0.30 mag (e.g. Hamuy et al. 1996), but this dispersion can be significantly 
reduced by using an empirical relationship between the light-curve shape and the raw ob- 
served peak magnitude. Phillips (1993) showed that the absolute magnitudes of SNe la are 
tightly correlated with the initial decline rate of the light-curve in the rest-frame S-band 
- intrinsically brighter SNe have wider, more slowly declining light-curves. Using a linear 
relationship to correct the observed SN peak magnitudes, parameterized via Amis {B), the 
decline in 5-band magnitudes 15 days after maximum light, the typical scatter in the rest- 
frame 5-band of local SNe can be reduced to ~ 0.17 mag (Hamuy et al. 1995, 1996; Phillips 
et al. 1999). An alternative technique, the Multi-color Light Curve Shape (MLCS; Riess 
et al. 1995, 1996, 2004a; Jha 2002), results in a similar smaller dispersion in the final peak 
magnitudes, particularly if multi-color light-curves are utilized in the fitting procedure. 

The approach used in this paper is to parameterize the SN light-curve time-scale via 
a simple stretch factor, s, which linearly stretches or contracts the time axis of a template 
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SN light curve around the time of maximum light to best-fit the observed light curve of the 
SN being fit (Perlmutter et al. 1997, 1999; Goldhaber et al. 2001). A SN's 'corrected' peak 
magnitude (m^ rr ) is then related to its raw peak magnitude (m^ w ) via 

m c ™ = + a (s - 1) (1) 

where s is the stretch of the light-curve time-scale and a relates s to the size of the magnitude 
correction. Sophisticated versions of this technique incorporating color corrections (e.g. Guy 
et al. 2005) lead to r.m.s. dispersions of ~ 0.15 mag in local SNe. 

2.2. SN model 

The apparent observed flux / of a SN la at redshift z through a filter with total system 
response T(A) is given by 

= tt~z l°° T{x)s (rb' *' s ' ET - y ' E *- y ) dA (2) 

where S(\,t, s, E™ v , E B °* V ), hereafter S, is the rest-frame flux distribution of the SN (in 
ergs -1 cm -2 A -1 ), which depends on the phase of the SN relative to maximum light (t, 
measured in rest-frame days), the SN stretch (s), the color excess due to the Milky Way 
along the line-of-sight to the SN (Eg™ v ), and the color excess of the SN host galaxy (E B ^y). 
S is normalized so that its intrinsic, unextincted, -B-band magnitude at maximum light 
(t = 0) m^ w is given by 

m -w = Mb + 5 \ogV c (z; n M , O0 - a(s - 1) (3) 

where V c = HqcIl is the "Hubble constant-free" luminosity distance, and Mb = M B — 
5 log H + 25 is the "Hubble constant-free" -B-band absolute peak magnitude of a SN la 
with s = 1 (e.g. Goobar & Perlmutter 1995; Perlmutter et al. 1999). In this model, we set 
a = 1.47 and M. B = —3.48 from the primary fit of Knop et al. (2003). (This formalism can 
of course be generalized to any SN type by removing the dependence on the stretch, s.) 

Model spectra and light-curves for the SNe are also required. The basic, normal SN la 
spectral templates we use are updated versions of those given in Nugent et al. (2002), which 
provide spectral coverage from 1000A to 25000A in wavelength, and from t = —19 to t = +70 
rest-frame days in the SN light-curve. The light-curve templates for calculating rest-frame 
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.B-band fluxes relative to maximum light (before any stretch-color correction) are the rest- 
frame light-curve templates of Goldhaber et al. (2001) and Knop et al. (2003), linearly 
interpolating between days as required. Light-curve information is provided for t — — 19 to 
t = +330 rest-frame days; for epochs before t = —19 days, the SN flux is assumed to be 
zero. 

We also incorporate the observed relationship between stretch and color from Knop 
et al. (2003), which describes the variation in SN rest-frame U — B, B — V, V — R and 
R — I colors as a function of stretch. For a given epoch, the spectral template is adjusted to 
have the colors for a given stretch using a spline which smoothly scales the spectrum as a 
function of wavelength. Extinction is applied to the altered template using the reddening law 
of Cardelli et al. (1989) updated in the near-UV by O'Donnell (1994), with R v = 3.1, first 
on the rest-frame spectrum with E^ty, and then on the redshifted spectrum with EgZy- 
We use the filter responses of Bessell (1990) and the a-Lyrae spectra energy distribution of 
Bohlin & Gilliland (2004) to perform a conversion to the AB system when required. 



For each SN candidate, we need to perform a fit to the n observed flux values O with 
uncertainty O err which depend on the filter T used for the observation, and the time of the 
observation U. We calculate how well a given model matches the observed data using the x 2 
statistic, i.e. 



where fx is defined in equation (2). Our task is to derive the model parameters z, s, t max , 
Eb°-v which best fit the data, i.e. which minimize \ 2 ■ While a grid search in x 2 -space is 
possible, in practice, with the number of free parameters in the model and typically a large 
number of potential SNe to fit, such a procedure becomes too time consuming for "real-time" 
operation. 

Instead, the best-fitting model templates are located using MPFIT 1 , a robust non- linear 
least-squares fitting program written in the IDL programming language translated from the 
MINPACK-1 FORTRAN package (Jorge et al. 1980). E%Z V is taken from the dust maps of 



http://cow.physics.wisc.edu/~craigm/idl/idl.html 



2.3. Fitting technique 
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Schlegel et al. (1998) and held fixed in the model. The output from this procedure are the 
best-fitting parameters of redshift, stretch and epoch, and optionally E^y. An additional 
parameter, termed a "dispersion magnitude" (dm) which is added to or subtracted from 
the effective 5-band magnitude of a model SN, is also fit. Physically, this accounts for 
the dispersion in light-curve-shape corrected peak magnitudes seen in local SN samples, of 
~ 0.17 mag; in our fits, this value is allowed to vary from -0.2 to 0.2 mag, but is not used 
further in the analysis. The other fit parameters are allowed to vary over the following 
ranges: 0.6 < s < 1.4, 0.05 < z < 1.4, and -0.10 < E^ v < 0.30, representing the typical 
parameter range of SNe la found in rolling searches. 

Once a first fit has been performed, the observed data are trimmed to only include data 
lying in the rest-frame to +30 days, the epoch to which stretch is most effective (Goldhaber 
et al. 2001). Once a first redshift estimate is made, data in an observed filter which cor- 
responds to rest-frame R or redder (where the current stretch technique is known to have 
limitations) are excluded from future iterations. For the Sloan Digital Sky Survey (SDSS) 
filter system, this implies that if the predicted redshift is less than 0.50, z' data are excluded 
from the second iteration fit. 



3. Testing on simulated survey data 

We now test the performance of the SN selection technique described above on simu- 
lated datasets, as a consistency check that our technique can successfully recover the input 
parameters of fake SNe la and determine their properties without any bias, and that a basic 
discrimination between SNe la and core-collapse SNe can be made. We generate both SNe la 
and core-collapse SNe in this simulation according to the design of the SNLS. 

3.1. SNLS survey design 

The Canada-France-Hawaii Telescope Legacy Survey (CFHT-LS) began in June 2003 
using the queue-scheduled wide-field imager "Megacam" on the CFHT on Mauna Kea. Mega- 
cam consists of 36 2048x4612 pixel CCDs covering a l°xl° field-of-view. CFHT-LS com- 
prises three separate surveys (Deep, Wide, and Very Wide), with the SNLS exploiting data 
from the Deep component. This component of the survey conducts repeat observations of 4 
low-extinction fields distributed in right ascension (Table 1). 

In a typical month, each available field is imaged on five epochs in a combination of 
g'r'i'z' filters which closely match those used for the photometric calibration of the SDSS. 
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r'i' observations are always taken, and g'z' observations are arranged according to the lunar 
phase. The typical break-down in exposure times for each field is given in Table 2; in 
general, each of the epochs are spaced 4-5 days apart (~3 days in a typical SN rest-frame). 
(u' data of the fields are also taken, but are not time sequenced as they are not used for 
the SN light-curves.) Each field is typically observed for 4-5 continuous months, plus a final 
month (with fewer epochs) to complete light-curves, resulting in around 20 "field-months" 
in a calendar year. The queue-scheduled nature of the observations provides a powerful 
protection against inclement weather, with the result that very few epochs are totally lost - 
most months have received the full allocation of 5 observation epochs per field. Consequently 
high-quality, multi-filter, continuous light-curves for each SN candidate lasting many months 
are routinely obtained. 



3.2. Simulating the SNLS 

We generate synthetic datasets as follows. For SNe la, events are drawn randomly from 
a population with truncated Gaussian distributions in redshift (z = 0.6, a = 0.4, z min = 0.15, 
z max = 1.05), stretch (s — 1, o — 0.25, s min = 0.6, s max = 1.4), dispersion in peak magnitude 
(dm = 0, o = 0.17), and host galaxy extinction {E%% = 0, a = 0.2, E\°% > -0.1). 

In contrast to SNe la, core-collapse SNe are quite heterogeneous events, with a large 
variation in their light-curves, spectra, colors and peak magnitudes; currently there is no well- 
established unifying scheme (c.f. stretch or Amis (B) for SN la) for the spectral and light- 
curve observations (though hints of a relationship between light-curve shape and luminosity 
have been observed for SNe Ic; Stanek et al. 2005). Nonetheless, we attempt to simulate 
them here using what little data is currently available in the literature, with the caveat that 
our simulations may not represent the full range of core-collapse SNe that could be observed. 

We generate four sub-types of core-collapse SNe, all based on spectral templates con- 
structed by one of us (PEN). For SN Ib/c, we use the spectral templates 2 described in 
Levan et al. (2005), and our spectral templates for SN IIL and SN IIP are similar to those 
used in Gilliland et al. (1999). We also created a SN Iln template based on SN 1999el (Di 
Carlo et al. 2002). The peak absolute magnitudes and dispersions for the core-collapse SNe 
are taken from Richardson et al. (2002); Gaussian distributions are again assumed. The 
other parameters defining the core collapse population are truncated Gaussians: redshift 
(z = 0.5, a = 0.3, z min = 0.05, z nmx = 0.6), and host galaxy extinction (-E^-V = 0, cr = 0.3, 
E^b-v > —0.2). These values differ from those used for the SNe la, and reflect their intrinsic 



2 The templates can be obtained at http://supernova.lbl.gov/~nugent/nugent_templates.html 



- 9- 



faintness when compared to SNe la, and likely increased dust extinction on these events (due 
to their association with recent star-formation). 

The simulations use an observational cadence similar to that of the SNLS i.e. 5 observa- 
tion epochs spread over 16 day-long observing periods centered around new moon (Table 2). 
Any given epoch has a 15% chance of being totally lost due to external factors, such as long 
periods of adverse weather. For each of these epochs, the rest-frame template SN spectrum 
is adjusted for the stretch-color relationship (for SNe la only), and additional dispersions in 
SN color are added, accounting for the natural variation in SN properties. For SNe la, the 
sigma of the dispersions used are 0.08 mag for U — B, and 0.03 mag for B — V, V — R and 
R — I. These dispersions are those intrinsic to SNe la after stretch-color relations have been 
applied. The increased dispersion in U — B over the other colors reflects the increasing sen- 
sitivity of the bluer passbands to chemical abundances (e.g. Hoflich et al. 1998, 2000; Lentz 
et al. 2000), and represent the broad range of observed values found in the literature (Nobili 
et al. 2003; Knop et al. 2003; Jha et al. 2005). For core-collapse SNe, the dispersion we use 
are: 0.20 mag for U — B, and 0.15 mag for B — V , V — R and R — I. The increased values are 
deliberately conservative, and represent our relatively poor knowledge of core-collapse SNe. 

The SN spectrum is then scaled to the expected effective £?-band magnitude for that 
phase (assuming an Qm = 0.30, Q\ = 0.70 cosmology), the stretch-luminosity relationship 
and the dispersion in peak magnitude from dm are applied, host extinction is applied, the 
spectrum is redshifted, and finally Milky Way extinction is applied. The spectrum is then 
integrated through the SNLS filters to generate an apparent magnitude in each filter, and the 
Megacam zeropoints used to calculate the counts in that filter given the exposure times listed 
in Table 2. The Megacam filter responses are those measured and provided by CFHT, and 
are additionally multiplied by the reflectivity of the CFHT primary mirror, the throughput 
of the wide-field corrector plus optics, and the average Megacam CCD quantum efficiency 
These effective filter responses can be seen in Figure 1, together with a typical SN la spectrum 
at a range of redshifts typical of SNLS candidates. 

Poisson noise, a filter-dependent sky noise (measured from real data in dark and gray 
sky conditions), and a filter-dependent systematic "flat-field error" noise is added to the 
predicted counts, and the counts converted back into fluxes using the original zeropoint, but 
with an additional zeropoint uncertainty (a = 0.05 mag) applied, typical of that present in 
real-time SNLS data (the zeropoint uncertainty for the final photometry is 0.02 mag or less, 
see Astier et al. 2005). Additional systematic noise is added to the z' (and a smaller amount 
in i') data to account for the imperfect fringe correction in real-time data. The result is a set 
of observations which should closely match real SNLS data. These synthetic SN observations 
are then passed to the SN selection code, and the code is run to locate the best-fitting SN 
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parameters. This experiment was repeated for 2000 synthetic SNe. The final population 
of SNe has the following percentage composition: SNe la 50%, SN Ib/c 25%, SN IIL 8%, 
SN IIP 8%, and SN Iln 8%. 

3.3. Simulation Results for SNe la 

We now examine the results of the fits to the simulated SN la photometry (the core 
collapse SNe are discussed in Section 5). In Figures 2 and 3 we compare the results of the 
input and recovered SN la parameters. This figure shows the comparison of actual versus 
recovered redshift, and the difference in the recovered stretch and dm parameters. We also 
show how the accuracy of the recovered parameters depends on the number of data epochs 
in the SN light-curve which are available at the time of fitting (i.e. 2, 3 or 4 pre-maximum 
light epochs). The mean and standard deviations of these distributions are shown in Table 3. 
Fits made to the entire light-curve were also run (i.e. including post-maximum light data); 
the results were very similar to Figure 2, with the exception that the distributions in stretch 
and dm are significantly tighter (see Table 3), as would be expected when there is more 
information available to constrain these parameters. 

For these synthetic, normal SNe la, the SN selection code appears capable of recovering 
their parameters without any apparent bias, even when only considering pre-maximum light 
data and when unknown zeropoint uncertainties and dispersions in the SN properties are 
included. Furthermore, the estimates of redshift are good to < 0.1 even when only two 
pre-maximum light-curve epochs are used in the fits. 

However, these SN la simulations do not address three important questions. First, do 
our synthetic SNe la successfully represent normal observed SNe la? Is the code capable 
of recovering all types of SNe la without any bias against certain parts of SN la parameter 
space? Though our simulations generate SNe la that should encompass the full observed 
range of low-redshift SN la properties, the possibility that some characteristics of SNe la are 
not represented, or that the code may reject them as SNe la, cannot be ignored. Secondly, is 
our selection technique capable of distinguishing non-SN la events (e.g., core-collapse SNe) 
from true SN la events on early-time data (of great importance in a rolling-search)? Thirdly, 
does the choice of the cosmology used in the light-curve fit (equation (3)) bias the SNe la that 
are selected for spectroscopic follow-up? We address these questions in the next two sections, 
using both the simulations described above, and using empirical data on spectroscopically 
confirmed SN events discovered via the SNLS. In Section 4, we describe the relevant aspects 
of the SNLS, and discuss the application of the SN selection code to the SNLS in Section 5. 
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4. Application to the SNLS 

In this section, we describe the application of the SN selection code to real-time data 
obtained as part of the SNLS. The results of the analysis, and a comparison with the simu- 
lations of Section 3, are then described in Section 5. 

4.1. SNLS Real-time data reduction 

A real-time reduction of the data are automatically performed by the CFHT-developed 
Elixir data reduction system (Magnier & Cuillandre 2004) 3 , which performs basic de-trending 
of the images, including bias-subtraction, flat-fielding and a basic fringe subtraction in the 
i' and z' filters. This real-time analysis is somewhat different to that of the final Elixir 
reductions released through the Canadian Astronomy Data Centre (CADC) 4 . Final processed 
data are able to use all data taken during a given queue run to construct master flat-fields and 
fringe frames, which is not possible in real-time where library versions of these calibration 
frames from the previous queue run are used. In particular, flat-fielding and, in the case of 
z' data, fringe subtraction are both inferior in real-time data. 

The Elixir-reduced data are then processed through two independent search pipelines, 
written by members of the SNLS collaboration in Canada and France, from which combined 
candidate lists are generated. The two pipelines typically produce candidate lists which 
agree at the 90% level down to %' ~ 24.0. In this paper, we refer exclusively to the Cana- 
dian pipeline, though the details of the French pipeline are similar in most respects. Full 
details of these reductions will be presented in forth-coming papers; here we present a brief 
overview. We perform a photometric alignment to secondary standards within each field 
using a multiplicative scaling factor, and astrometrically align the data in different filters 
using custom-built astrometric reference catalogs. Finally, the frames are re-sampled to a 
common pixel coordinate system, and a median stack in each filter is generated, rejecting 
cosmic-ray events and other chip defects. 

The point-spread functions (PSFs) of the image stack and the deep reference image 
constructed from observations from the previous year are matched using a variable kernel 
technique, and a subtraction image generated containing only sources which have varied since 
the reference epoch. This subtraction image is searched for SN candidates using automated 
techniques, with all likely candidates screened by a human eye to ensure obviously non-stellar 



3 http : //www . cf lit . hawaii . edu/Instruments/Elixir/ 
4 http : //cadcwww . dao . nrc . ca/ 



- 12 - 



sources are not included on the potential candidate follow-up list. 

PSF-fit and aperture photometry for all new and previously detected sources is measured 
and the flux information is stored in a publicly- available database of all SN candidates 5 . New 
candidates are automatically cross-correlated with the database of existing variable sources 
(to identify previously known AGN and variable stars), and back-tracked through previous 
epochs in order to construct as much information as possible as to their nature. Additionally, 
each candidate is given a provisional category of either likely supernovae (SN), likely AGN 
(if the detection is precisely centered and appears intermittently over long periods of time), 
likely variable star (if the candidate host appears stellar and PSF-like), or likely moving 
object (if the candidate moves between epochs). Where any reasonable doubt exists, the 
SN classification is retained. The result is a database of flux information for every variable 
candidate detected in the SNLS. 

Spectroscopic follow-up for the SNLS is currently provided by dedicated follow-up pro- 
grams on the European Southern Observatory Very Large Telescope (VLT; 60 hours/period), 
Gemini North and South (60 hours/semester) and Keck-I/II telescopes (4 nights in the "A" 
semester of each year), as well as other programs at Keck-I making detailed studies of select 
SNe (4-8 nights a year). During the time-frame covered by this paper, the VLT observed 61 
candidates (Basa et al., 2005 in prep.), Gemini observed 40 candidates (Howell et al. 2005) 
and Keck observed 21 candidates, 7 as part of routine follow-up, and 14 for a detailed study 
program (Ellis et al., in prep.) Some candidates were observed on more than one occasion. 
Full details of the reductions and SN typing methods employed can be found in the papers 
referenced above. Typically, redshifts are determined from host galaxy spectral emission or 
absorption features, or via a x 2 template matching technique to a large library of local SN 
spectra of all types (Howell et al. 2005). 



4.2. Analysis 

The photometry and candidates that we use in this analysis are taken from the Canadian 
real-time analysis corresponding to the period May 2004 to January 2005 - or approximately 
3 "field-seasons" . Over this period, the SN selection code was used to assist in the selection 
of spectroscopic follow-up candidates. 440 potential SN candidates were located by our real- 
time analysis software (excluding likely AGN and variable stars), of which 121 were followed 
spectroscopically. Of these, 70 have since been identified as SNe la or probable SNe la 
(denoted SN la* hereafter), 6 as being consistent with an SN I, 11 as SNe II, 4 as SN Ib/c. 



5 http : //legacy . astro . utoronto . ca/ 
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30 candidates remain unidentified or are non-SN events. During this period, in order to test 
the SN selection code, in some months we deliberately followed events that were predicted 
not to be SNe la to ensure we were not biasing against particular classes of SN la events. 

For this analysis, we run the SN selection code on all SN candidates discovered during 
this period (fixing Eg w v as appropriate for each candidate), firstly fitting on the entire 
real-time light-curve, and secondly fitting just the data that were available at the time 
the follow-up decision was made. This exact date was not accurately recorded or easily 
reconstructed for all events; we approximate this by cutting-off the real-time light-curves 
at 2 days before maximum light in the SN rest-frame. Clearly, this cut restricts us to only 
consider candidates with at least 2 epochs prior to 2 days before maximum light, ~ 95% of 
the total SN la sample; furthermore, the approximation may lead to some fits using more 
data (and some less) than was actually used at the time the follow-up decision was made. 

5. Results 

5.1. Discriminating core-collapse SNe from SNe la 

We first examine the non-SN la events located by the SNLS. One important requirement 
of the SN selection technique is to not only return the most likely parameters describing a 
SN based on the fitting the real-time light-curve to a normal SN la model, but also to return 
some guide as to the likelihood that the target is indeed a normal SN la. As fits to non-SNe la 
should be of lower quality, we examine the x 2 °f the light-curve fits to examine whether it 
can be used as a reliable indicator. 

The most common form of contamination in our sample comes from core-collapse SNe. 
While AGN and variable stars are obviously present, they are usually easily screened out as 
they show significant variation on the time-scale of a year, and are usually already present 
in our database of variable objects. Core collapse SNe - the various sub-classes of SNe II 
and SNe Ib/c - are heterogeneous; the range in peak magnitudes (Richardson et al. 2002), 
spectra (Matheson et al. 2001), and colors can be considerable. The general trend however 
is that for a core-collapse SN and a SN la located at the same redshift, the core collapse SN 
will be fainter than the SN la. Typical SNe II are 1-2.5 mag fainter at peak than SNe la, 
and typical SNe Ib/c around 1.5 mag fainter (see table 1 in Richardson et al. 2002), though 
there are examples of more luminous or peculiar events in the literature (e.g. Germany et al. 
2000; Turatto et al. 2000). 

Our SN la selection technique uses both this peak magnitude information and the color 
of the observed light-curve, fitting a given candidate for redshift, phase, stretch, extinction, 
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and dispersion in the peak magnitude under the assumption that the object is a SN la in 
an assumed cosmology. This places a constraint on the redshift that can be fit (as SNe la 
are calibrated candles and have similar intrinsic luminosities), and hence restricts the range 
of observed colors that are consistent with a SN la as a function of SN brightness. We 
demonstrate this in Fig. 4, where we show the redshift evolution in the color of a SN la and 
two example models for a SN II (in this case an SN IIP and an SN Iln), and Fig. 5, where we 
show the evolution of SN colors as a function of phase. The SN la template is that discussed 
in Section 2.2 and the core-collapse templates are described in Section 3.2. 

Consider a SN IIP at maximum light located at z — 0.25 with g' — r' ~ and r' ~ 23. 
Though the color is also consistent with a maximum-light SN la at the same redshift (Fig. 4), 
the peak magnitude of the SN IIP is around 2 magnitudes fainter than a SN la would be 
(r' ~ 21). This r' ~ 23 peak magnitude is consistent with a SN la at z ~ 0.55; however, at 
this redshift the g' — r' color of a normal SN la is ~1 mag redder. Hence, SNe II (and other 
core-collapse events) can potentially be screened out by searching for objects that appear too 
blue in g' — r' for the fitted redshift when compared to the SN la model, or more precisely, 
rejecting objects that have a large x 2 P er degree of freedom (d.o.f) in the (/-band from the 
SN la fit. A measurement of color at exactly maximum light is not required; the same trend 
is apparent up to about 20 days after explosion (Fig. 5). 

The distributions of x 2 in the g'-band from our pre-maximum light fits to the simulated 
and real SNLS candidates are shown in Fig. 6. Example light-curve fits to both real SNe la 
and real core collapse SNe are shown in Fig. 7. As expected, the SNe la and core collapse SNe 
have quite different distributions, both for the simulated survey data, and for the real survey 
data. SNe la typically have a small x 2 ; the median value is 0.85 (1.01 for the simulated 
data), and 90% of SNe la have x 2 < 2.80 (x 2 < 2.91 in the simulations). Only 1% (0.9% in 
simulated data) have x 2 > 5. 

The x 2 f° r core collapse SNe is typically larger, with 10 of 11 events having a x 2 > 
5; all SNe II have x 2 > 5. The numbers for the fits to the simulated SNe are similar; 
80% of simulated core-collapse SNe have a x 2 > 5. Thus, the x 2 in the observed g' filter 
seems an efficient discriminant with which to remove the majority of core collapse SNe from 
spectroscopic follow-up samples. 

5.2. Redshift and phase precision 

The accuracy to which redshift and phase can be estimated are important quantities; 
the first allows an optimal spectrograph setting to be used, and the second ensures that SNe 
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are observed at a phase when the different sub-types show the most diversity and with the 
greatest contrast over their host galaxy. Here we investigate how well the predicted redshift 
and phase of the SNe la agree with actual values. 

The photometric-redshift (%,hot) versus spectroscopic-redshift (-2 sp cc) comparison for all 
spectroscopically observed objects where a z spcc could be determined is shown in Fig. 8. The 
plot is based on fits to just the data that was available at the time the follow-up decision 
was made (taken as pre-maximum light data). 

The agreement between z p hot and z spcc for SNe la is remarkably good. The median 
\Az\ /(l + z S pc C ) is 0.031, with 90% of the SNe la having |Az| /(l + z spcc ) < 0.08. The 
agreement for other SN types is, as expected, much poorer, as the differing colors, light- 
curves shapes and brightnesses of these SN leads to an incorrect redshift estimate. However, 
most of these core collapse candidates are rejected by the code as being SNe la due to 
the x 2 °f the fits; such rejected SNe are marked on the figure. As a comparison, Fig. 9 
shows the z p hot versus z spcc for fits based on the entire SN light-curve; here the median 
|Az|/(l + %,cc) is 0.025. Note that the apparent high accuracy of this technique does 
not mean that spectroscopic redshifts are not needed for cosmological analysis; clearly as 
a cosmology is assumed in order to derive the fitted redshift, this redshift could not then 
be used to re-derive cosmological constraints. Furthermore, as the nature of the events are 
not known to be SNe la a priori, and the assumed cosmology is a very useful prior when 
identifying the core collapse SNe, the technique could not be used to derive cosmological 
constraints (as can be done in some low-redshift SN la samples, e.g. Barris & Tonry 2004) 
without a significant improvement in our knowledge of the color evolution and light-curve 
morphology of core collapse SNe. 

The distribution of the predicted date of maximum light versus the date of maximum 
light based on the entire light-curve is shown in Fig. 10 (as before, predicted values are 
only based on candidate light-curves up to 2 rest-frame days before maximum light). The 
agreement is good, with 50% of the predicted times of maximum light lying within ±1.75 
rest-frame days of the actual date of maximum light, and 90% within ± ~ 4 days. 

5.3. Effect of assumed cosmology 

In this section, we investigate any possible bias in the population of objects that would 
be selected for spectroscopic follow-up based on the "assumed" cosmology that is used in 
the light-curve fit. We repeat the SN selection fits on the real survey data using a different 
cosmological model, and compare the identity of the rejected candidates with those rejected 
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when using the standard cosmology. 

Subtle changes in the cosmology, for example varying the value of the equation of state 
parameter w, have no effect on the nature of the candidates that are rejected. For example, 
the candidates that are rejected (or selected) with w — — 1 are also rejected (or selected) 
with w = 0.8 or w = —1.2, as the differences in the apparent-magnitude/redshift relation 
are very small. A more drastic test is to compare to a cosmology that makes very different 
predictions in the apparent magnitude of an object as a function of redshift. We compare to 
an Einstein-de-Sitter (EdS) cosmology (Qm — 1-0, — 0.0); the apparent magnitude of a 
SN la differs by ~ 0.5 mag at z — 0.6. Any bias in the objects that are selected for follow-up 
by assuming an accelerating universe cosmology should be evident after these fits. 

In Fig. 11 we show the same -2 p hot versus z spcc comparison as in Fig. 8, but for the EdS 
Universe. There is a clear offset apparent in these plots; in an EdS Universe, the z p hot is 
found to be systematically larger than the -2 spcc , and the size of this effect increases at higher 
redshift. This is simple to understand; in a EdS universe, objects appear brighter than in 
a A-cosmology at the same redshift. If we live in a A-dominated Universe, the SN selection 
code will over-estimate the redshift (i.e. will try to make the objects fainter) if it assumes 
an EdS cosmology. 

The objects that are rejected based on the g' x 2 are the same when assuming an EdS 
Universe as when assuming a A-cosmology i.e. the assumed cosmology does not place a strong 
constraint on the objects that are selected for spectroscopic follow-up. However, assuming 
a A-cosmology does improve our estimates of candidate redshift s and hence stretches and 
observer-frame ages, which is essential to ensure that candidates are observed at the optimal 
time (Howell et al. 2005). The assumed cosmology does not provide the key discriminant 
in rejecting non-SN la events, but does enable a more efficient follow-up observation to be 
performed. 

5.4. Implications for cosmological measurements 

We now examine the SNe la that are rejected in any cosmology, i.e. SNe la that 
have sufficiently discrepant stretches, luminosities, or colors when compared to the standard 
template that is used in the fits, or SNe where the selection code simply fails to find the 
correct solution. If similar sub-types of SNe la were consistently rejected, this could introduce 
a bias in the types of SNe la that were sent for follow-up, with a possible consequent bias 
on the cosmological results that would be derived from the resulting spectroscopic sample. 
Alternatively, if the occasional rejections of a SN la were random, and related to the phasing 
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of the observations or the poor quality of a given epoch of real-time data, then the impact 
on the derived cosmological parameters would be very small. 

We investigate this effect using a simulated dataset similar to that in Section 3. We 
generate 700 high-redshift (0.15 < z < 0.9) SNe la (the SNLS end-of-survey goal), and 
combine with 300 local SNe la (0.01 < z < 0.12), a sample size equivalent to that likely to 
be available at the end of the SNLS from projects such as the Carnegie Supernova Project 
(Freedman 2005) and the Nearby Supernova Factory (Aldering et al. 2002). The 700 high- 
redshift SNe are fit using the SN selection code, the pre-maximum light data fit, and those 
SNe which are formally rejected are noted. Typically 1-2% are formally rejected by the code, 
usually because of a lack of data before maximum light (due to weather in our simulations) 
rather than an intrinsic property of the SN. We then perform a full cosmological fit of w 
and Ojtf, firstly on the full 700+300 SN la sample, and then on the same sample less the few 
high-redshift SNe that were rejected by the SN selection code based on the pre-maximum 
light fits. The cosmological fits follow a similar method as that used in Astier et al. (2005), 
and produce confidence contours in w and Qm- In all cases, the input cosmological values 
were Qm = 0.3 and w — — 1, and a flat Universe was assumed. 

As each survey simulation generates different SNe, we repeat the full simulation four 
times; the confidence contours for the recovered cosmological parameters in Qm and w are 
shown in Fig. 12. Constraints from the baryon acoustic oscillations (BAO) are also included 
(Eisenstein et al. 2005). The differences between the cosmological fits for the different sam- 
ples are very small. In all cases, the difference in the derived value of w is less than 0.001 
(~0.02<r) when including the BAO constraints, and less than 0.015 (~0.05a) when the BAO 
constraints are not used. We therefore conclude that the SN selection method adds no 
significant bias to the eventual cosmological results of the SNLS. 

5.5. Potential SN la targets not followed spectroscopically 

An alternative, more empirical, method of how successful the SN selection code is at 
locating candidates for follow-up in real-time, is to examine the fraction of good SN la candi- 
dates rejected by the code. The impact on the determination of the cosmological parameters 
of these SNe was shown to be negligible in the previous section; however, this is also im- 
portant for future studies of the SN la that rely on near-complete spectroscopic samples of 
SNe la being available (e.g., SN rates). We assess this by examining the distribution of the 
fitted SN stretches of our candidates once their light-curves are complete. Virtually all nor- 
mal SN la stretches lie in the range 0.7 to 1.3 (e.g. Perlmutter et al. 1999; Astier et al. 2005). 
SNe la with very different stretches have been discovered locally, for example the broad 
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s = 1.6 light-curves of very rare SN 2001ay-like SNe (Howell & Nugent 2004) or the faint, 
fast declining SN 1991bg-like objects (e.g. Filippenko et al. 1992; Garnavich et al. 2004). Al- 
though SNe such as SN 1991bg have been shown to be calibrateable as standardized candles, 
and have been used in determinations of the Hubble constant (Garnavich et al. 2004), these 
classes of peculiar SNe are unlikely to be useful for cosmological studies at high-redshift - 
they may not follow the lightcurve-shape/luminosity relationship (in the case of SN 2001ay), 
and may be naturally selected against at high redshift due to their faintness (in the case of 
SN 1991bg-like SNe, which are around 2 magnitudes under-luminous). In "distance-limited" 
surveys, which likely provide a census of all types of SNe la, around 65% of SNe la appear 
spectroscopically "normal" (Li et al. 2001), a number which is probably a lower limit on 
the number of cosmologically useful SNe, as it excludes SN 1991T-like events that may still 
be calibrateable candles. This fraction will only rise in high-redshift searches as the fainter 
SN 1991bg-like SNe are selected against. Hence candidates with fitted stretches outside of 
the range 0.7 to 1.3 are either unlikely to be SNe la, or unlikely to be cosmologically useful 
SNe la, while those with 0.7 < s < 1.3 are likely to be useful SNe la. 

Figure 13 shows the distribution of fitted stretch for all candidates (even if the spectro- 
scopic redshift is known, redshift is not fixed in these fits). As expected, the stretches for 
the spectroscopically confirmed SNe la lie in the expected range 0.7 to 1.3. Furthermore, of 
the candidates not selected by the code for spectroscopic follow-up, few appear to be good 
candidates for normal SNe la once objects for which a spectroscopic observation would be 
too challenging have been removed. Objects which were not possible to observe spectroscop- 
ically are defined as those for which the peak %' magnitude is > 24.25, or for which the peak 
i! percentage increase over the host galaxy is < 30%. (SNe that do not meet these cuts are 
rarely positively identified - see Howell et al. 2005) Fig. 14 shows the same objects broken- 
down into the same categories, but the histogram is of the fitted redshift rather than the 
stretch. As would be expected, objects defined as too difficult to attempt spectroscopically 
are only found at z > 0.65, where the contrast over the host galaxy is typically smaller as 
the hosts have a smaller apparent size on the sky - the angular separation between SN and 
host is smaller, and the SN light becomes more mixed with the galaxy light. At z < 0.65, of 
all the total SN candidates that were consistent with a SN la, only ~ 5% were not followed 
spectroscopically where an observation would have been possible. 

6. Summary 

In this paper, we have presented a selection technique for high-redshift supernovae 
searches that can be used to identify SNe la after only 2-3 epochs of multi-band photometry. 
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Using both simulated SN data and the SNLS real-time data as a test-case, we have shown 
that the technique is able to discriminate between SNe la and core collapse SNe (SNe II, 
SNe Ib/c) based on the quality of the fit in the g' filter. For SNe la, the technique accurately 
predicts the redshift, phase and light-curve parameterization of these events to a precision 
of |Az| /(l + <2 S pcc) < 0.09 in redshift, and ±2-3 rest-frame days in phase, and there is no 
apparent bias on cosmological parameters derived using SNe la selected in using this method. 

The technique is now routinely used within the SNLS to help select priority candidates 
for spectroscopic follow-up and confirmation as SN la for use as standard candles in the 
derivation of cosmological parameters (Astier et al. 2005); the improvement in spectroscopic 
success the method brings is discussed elsewhere (Howell et al. 2005). These techniques will 
be used in other current and future planned surveys, which are likely to be "spectroscopically 
starved", with many more candidates than can be followed up. 

The obvious extension to this method will be when "final" data (versus the real-time 
data used here) becomes available for SNLS. With many of the real-time uncertainties re- 
moved, and a high quality fringe subtraction capable of removing systematics, this SN selec- 
tion technique will provide a logical method for determining the types and redshifts for the 
many hundreds of SNe that it is not possible to observe spectroscopically, particularly when 
combined with photometric redshifts for the host galaxies that will be routinely available 
for the CFHT-LS Deep Fields. Such a combination of data would allow a calculation of the 
rates of all SN types out to z = 1. 
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were also indispensable in making possible real-time data reduction at CFHT. Based on ob- 
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(NRC) of Canada, the Institut National des Science de l'Univers of the Centre National 
de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work 
is based in part on data products produced at the Canadian Astronomy Data Centre as 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative project of NRC 
and CNRS. Canadian collaboration members acknowledge support from NSERC and CIAR; 
French collaboration members from CNRS/IN2P3, CNRS/INSU and CEA. Based on obser- 
vations (Program-IDs: GN-2004B-Q-16, GS-2004B-Q-31, GN-2004A-Q-19, GS-2004A-Q-11, 
GN-2003B-Q-9, and GS-2003B-Q-8) obtained at the Gemini Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under a cooperative agree- 
ment with the NSF on behalf of the Gemini partnership: the National Science Foundation 
(United States), the Particle Physics and Astronomy Research Council (United Kingdom), 
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the data presented herein were obtained at the W.M. Keck Observatory, which is operated 
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Table I. CFHT-LS Deep field locations 



Field 


RA (J2000) 


DEC (J2000) 


Other data 


Dl 


02:26:00.00 


-04:30:00.0 


XMM-Deep, VIMOS, SWIRE, GALEX, VLA 


D2 


10:00:28.60 


+02:12:21.0 


COSMOS/ACS, VIMOS, SIRTF, GALEX, VLA 


D3 


14:19:28.01 


+52:40:41.0 


(Groth strip); DEEP-2, SIRTF, GALEX 


D4 


22:15:31.67 


-17:44:05.7 


XMM-Deep, GALEX 



Table 2. Exposure times and epochs for a typical month in SNLS 





Epoch 1 


Epoch 2 


Epoch 3 


Epoch 4 


Epoch 5 


~Days w.r.t. new moon: 


-8 


-4 





+4 


+8 


g' a 


5 x 225s 




5 x 225s 


5 x 225s 




T> 


5 x 300s 


5 x 300s 


5x300s 


5 x 300s 


5 x 300s 


%' 


7x520s 


5 x 360s 


7x520s 


5 x 360s 


7x520s 


z' 


10 x 360s 




10 x 360s 




10x360s 



a For the g' filter, epochs 1/2 and epochs 4/5 can be exchanged depending on moon position. 
From August 2005, 4 g' epochs have been obtained. 
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Table 3. "actual minus recovered" SN parameters for the 2000 synthetic SNe la 



Parameter 


Simple 


Outlier resistant 




Mean 


a 


Mean 


a 


Pre-maximum light- 


-curve 


epochs only: 




% ^recovered 


0.098 


0.240 


0.007 


0.025 


S ^recovered 


0.204 


0.423 


0.015 


0.131 


d.T/2 CLTTZ rccoverec [ 


-0.004 


0.144 


-0.004 


0.140 


Entire light-curve: 










% ^recovered 


0.006 


0.041 


0.007 


0.022 


S ^recovered 


0.041 


0.099 


0.021 


0.044 


d.771 d.TTZ recoverec [ 


0.010 


0.142 


0.012 


0.132 
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Fig. 1. — The mean effective responses of the g', r', i', and z' filters used in CFHT-LS, 
shown shaded left to right in the figure. The filter responses are those measured by CFHT, 
additionally multiplied by the reflectivity of the primary mirror, the throughput of the wide- 
field corrector optics, and the typical Megacam CCD quantum efficiency Plotted underneath 
is a maximum-light s = 1 SN la spectrum at redshifts of 0.3 (dotted), 0.6 (solid) and 0.9 
(dashed), spanning the range of SNe la followed by the SNLS. The SN template is a revised 
version from Nugent et al. (2002). 
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Fig. 2. — Results of the SN selection code run on a synthetic dataset of 2000 SNe la generated 
with observational parameters similar to that found in the SNLS (see § 3.2 for simulation 
details). This figure shows the comparison of actual and recovered redshift, and how this 
varies depending on the number of pre-maximum light-curve epochs are used in the fit, 
excluding pre-SN explosion epochs, from 2 epochs (top panel) to 4 epochs (lower panel). 
Only pre-maximum light data (including epochs of non detections) are used when recovering 
the redshift. On each epoch, 3 filters are observed, either g'r'i' or r'i'z'. 
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Fig. 3. — As Figure 2, but showing the distribution of actual minus recovered stretch (left 
panels), and for dm (right panels). 
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Fig. 4. — The observed g' — r' and r' — i' colors of SNe la compared to some SNe II models 
as a function of redshift. As SNe la typically appear 1-2 mag brighter than SNe II located 
at the same redshift, the strong evolution in SN la g' — r' color with redshift provides a 
strong discriminant against core-collapse SNe, which appear too blue for a redshift based on 
a SN la fit. 
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Fig. 5. — The observed g' — r' and r' — %' colors of SNe la compared to some SNe II models 
as a function of rest-frame phase of the SN. Two redshifts are shown; z = 0.25 and z = 0.55. 
Note that SNe II and SNe la peak at different times after explosion; these are indicated with 
different vertical lines in the figure. 
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Fig. 6. — The distribution of g' x 2 P er d.o.f. for the fits of all SN candidates for the 
simulations (Panel a, left; Section 3) and for real SNe discovered by SNLS (Panel b, right). 
In panel a) 'Core collapse SNe' represent all the various sub-types of core-collapse SNe as 
described in the text. In panel b), 'Confirmed SNe la' are those SNe identified as either 'la' 
or 'la*' from their spectra, 'Confirmed core collapse SNe' are those identified as any variant 
of a SN II, 'Other spectra' includes any other object observed spectroscopically regardless 
of the type of object. 
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Fig. 7. — Example light-curve fits to early real-time photometry. In all cases the light-curves 
are fit to a SN la model, with the phase, redshift and stretch free parameters in the fit as 
described in the text. The top row shows fits to known SNe la; from left to right 04D4in 
(z = 0.52), 04D4jr (z = 0.48), and 04Dlow (z = 0.92). The bottom row shows fits to 
known core-collapse SNe; from left to right: 04Dlln (SN IIP; z = 0.21), 04Dlme (SN IIP; 
z = 0.30), and 04Dlla (SN Ib/c; z — 0.32). In all cases, the deviations from the SN la model 
are significantly larger for the core-collapse SNe, allowing an efficient identification of such 
objects based on the g' x 2 - 
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Fig. 8. — The comparison between 2 p hot and z spcc for the candidates followed as part of the 
SNLS. The photometric redshifts are from the SN and are based on pre-maximum-light data 
only. The SN are coded according to the type of the event, and those events that would 
have been removed from our follow-up sample based purely on the g' x 2 (P er d.o.f.) from the 
SN la fit are boxed with a diamond. The two rejected SNe la have only poor quality pre- 
maximum light data; in reality more data was utilized for these candidates before a follow-up 
decision was made. The z p hot errors are those reported by the least-squares minimizer. The 
over-plotted dashed lines show a 1:1 correspondence between 2 p hot and z spcc as well as ±10% 
in (1 + z). The assumed cosmology is Vt M = 0.25, = 0.75, and the SNe la include both 
'la' or 'la*' events. 
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Fig. 9. — As Fig. 8, but using the entire light-curve in the fit rather than just the pre- 
maximum light data. Note that more SNe appear on this plot than on Fig. 8, as this plot 
includes SNe for which there is no pre-maximum light data. 
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Predicted phase - fit phase (SN rest-frame days) 

Fig. 10. — The distribution of the predicted time of maximum light minus the actual time of 
maximum light in the SN rest-frame. The "predicted time of maximum light" is based only 
on pre-maximum light light-curve data i.e. the data available at the time of spectroscopic 
follow-up. The "actual time of maximum light" is based on the entire SN light-curve. The 
vertical lines denote the region containing 50% (dashed) and 90% (dotted) of all candidates. 
Only confirmed SNe la (or SNe la*) with spectroscopic redshifts are shown. The uncertainty 
on the time of maximum light from fits based on the entire light-curve is ~ 0.5 — 1.0 days. 
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Fig. 11. — As Fig. 8, but assuming an Q M = 1.0, = 0.0 cosmology. Changing the 
assumed cosmology does not alter the SNe that are rejected as being SNe la by the fitting 
code, though the agreement between z spec and z p hot is poorer. 
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Fig. 12. — The cosmological fits in w and VLm to four simulated SNLS-like SN samples with 
and without the few SNe la formally rejected by the SN selection code. Each panel represents 
a different SNLS-like survey; the differing positions of the contours is expected and arises 
from random variation from simulation to simulation. In each panel, the solid contours 
indicate the full sample of 700 SNLS SNe plus 300 low-redshift SNe, the dotted lines are the 
same sample minus those high-redshift SNe rejected by the SN selection code. In each case, 
the confidence contours with and without the baryon acoustic oscillations constraint (dashed 
line; Eisenstein et al. 2005) are shown. The number in each panel indicates the number of 
high-redshift SNe not rejected by the selection code. 
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Fig. 13. — The histogram of fitted stretch for all SN candidates from the SNLS during the 
time period covered by this paper. 'Confirmed SNe la' are those SNe identified as either la or 
la* from their spectra, 'Other spectra' includes any other object observed spectroscopically 
regardless of the type of object that was not identified as a SN la, 'Rejected on x 2 of g' 
fit' includes all remaining objects that are rejected by the SN selection code based on the 
fit to the (/-band light-curve, 'Not suitable for spectroscopy' includes all remaining objects 
which at maximum light are either too faint or have too low percentage increase for a likely 
successful spectroscopic observation, and 'Remaining potential SNe' includes all remaining 
objects which are not in the above categories and which are potentially SNe la not observed 
spectroscopically. The photometric parameters are based on photometric-redshift fits to the 
entire candidate light-curve. The peaks at s = 0.6 and s = 1.4 reflect the range over which 
stretch is allowed to vary in the fits (0.6 < s < 1.4). 
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Redshift fitted to SN la model 

Fig. 14. — The histogram of SN la photometric redshifts for all candidates from the SNLS 
during the time period covered by this paper. The categories of objects are the same as in 
Fig. 13. The photometric redshift is based on fits to the entire candidate light-curve under 
the assumption that the object is a SN la; clearly for non-SNe la the z p hot will be inaccurate. 
The figure shows that there are few good candidates lying at %>hot < 0.65 which were rejected 
by the selection code in the real-time analysis. Note that the redshift distribution of the 
candidates 'Rejected on x 2 of g' fit' is not representative of the true redshift distribution of 
these events as these candidates are unlikely to be SNe la. 



